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Topological insulators in the presence of strong Coulomb interaction constitute novel phases of 
matter. Transitions between these phases can be driven by single-particle or many-body effects. On 
the basis of ab-initio calculations, we identify a concrete material, i.e. Ca 2 Pt 04 , that turns out to be 
a hole-doped weak topological insulator. Interestingly, the Pt -d orbitals in this material are relevant 
for the band inversion that gives rise to the topological phase. Therefore, Coulomb interaction 
should be of importance in Ca 2 Pt 04 . To study the influence of interactions on the weak topological 
insulating phase, we look at a toy model corresponding to a layer-stacked 3D version of the Bernevig- 
Hughes-Zhang model with local interactions. For small to intermediate interaction strength, we 
discover novel interaction-driven topological phase transitions between the weak topological insulator 
and two Dirac semimetal phases. The latter correspond to gapless topological phases. For strong 
interactions, the system eventually becomes a Mott insulator. 
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Two classes of materials, which are in the center 
of focus of condensed-matter physics, are topological 
insulators (TIs) and strongly correlated systems. So 
far, their studies have taken rather disjunct paths: the 
striking feature of the TIs, namely the existence of 
stable (Dirac) surface/edges states in a bulk insulat¬ 
ing system, has first been predicted on the basis of 
“non-interacting” band models [1-5]. In the case of 
the “Bernevig-Hughes-Zhang” (BHZ) model [6], it has 
then triggered the first experimental realization of a 
two-dimensional (2D) quantum-spin-Hall (QSH) system 
based on CdTe/HgTe/CdTe quantum wells [7]. In this 
model, as well as in other semiconductor materials [8-10] 
a “band-inversion” mechanism is crucial, where conduc¬ 
tion and valence bands are inverted by spin-orbit cou¬ 
pling (SOC). This “non-interacting” path was soon car¬ 
ried over to 3D TIs [11-16], where the search for the 
corresponding 2D massless Dirac states is additionally 
spurred by the quest for a deeper understanding of topo¬ 
logical states of matter as well as wider spintronics ap¬ 
plications. 

Notable exceptions in this path to TIs are a recent im¬ 
plementation of local correlation effects into the 2D BHZ 
model, the Kane-Mele-Hubbard model as well as propos¬ 
als for a Kondo TI state [17-27]. In the study of the 
2D BHZ model, some of us showed that a system with 
topologically trivial parameters in the absence of inter¬ 
actions, i.e. a band insulator, can be driven into a QSH 
phase by local electronic correlations [28]. Interestingly, 
this transition can even become of first-order [29]. 

Here, we demonstrate on the basis of density- 
functional theory (DFT) that a weak TI phase can be 
achieved in a concrete 5 d transition-metal compound, 
i.e. hole-doped Ca 2 Pt 04 . This discovery stimulates us 
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to study electronic correlation effects embedded in a toy 
model for 3D weak TIs, which is a many-body general¬ 
ization of the layered BHZ model. We expect that the 
concrete material example combined with our many-body 
study of the layered BHZ model will provide a rich ma¬ 
terial playground for a deeper understanding of the com¬ 
petition between topology and interactions, in particular, 
in weak TIs. Somewhat surprisingly, we predict a topo¬ 
logical phase transition from a weak TI phase to a Dirac 
semimetal (DSM) phase that is driven by interactions. 
The DSM phase is characterized by 3D bulk Dirac points 
that are of topological origin. The reason is that their 
existence (even in the non-interacting limit) is forced by 
the change of a 2D topological invariant as a function of 
one wave vector, e.g. k z , if the 2D invariant is calculated 
for a fixed value of k z . 

When the most general on-site electronic correlations 
are included in our model via a DMFT calculation, we 
obtain a rich physics for small to intermediate interac¬ 
tions: Both the weak TI and the DSM phases are stable 
and the transitions between them can even be driven by 
electronic interaction. To the best of our knowledge, this 
is the first prediction for transitions between a trivial 
band-insulator to a DSM phase to a weak TI triggered 
by Coulomb interactions. 

Material example: In the first part of this Letter, we 
provide an example for a weak topological insulating 
phase in a particularity exciting 3D system. When exam¬ 
ining the experimental possibility of realizing an interact¬ 
ing topological phase, electronic correlations and strong 
SOC should ideally coexist. Most of the TIs discovered in 
3D are, however, semiconductors, whose nontrivial topol¬ 
ogy is generated by the strong SOC of heavy elements, 
such as bismuth or mercury. On the other hand, strong 
electronic correlations usually appear in the incompletely 
filled d- or /-electron shells with narrow energy bands. 
Thus, it would be desirable to find a nontrivial topology 
from a d — p band inversion. Then, electronic correla- 
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tions and topology could both matter. In fact, such a 
transition can appear in 5 d transition-metal-oxides. In 
this respect, perovskite iridates [30-33] have been theo¬ 
retically predicted to be a possible host for a correlated 
TI phase. Here, we propose that platinum oxides can be 
another promising platform for studying correlated TIs. 
Compared to iridates, this is rather unexplored territory. 
The example discussed below, i.e. the hole-doped tran¬ 
sition metal compound Ca2Pt04, shows that a weak TI 
phase can appear in platinum oxides which has not been 
found experimentally in iridates so far. 



FIG. 1. (Color online) Electronic structure of Ca 2 Pt 04 along 
the high-symmetry direction N(|,0,0)-r(0,0,0)-N(0,|,0)- 
X(i,i,0)-Z(i,i,i)-N(|,0,|)-X(0,0,|)-N(0,i,i). The light- 
yellow region shows the topological band gap which contains 
two inversions at T and Z between the bands with mainly 
0-p z (red circle) and Pt -d xy (green circle) characters. 

In Fig. 1, we display our electronic structure results 
of Ca2Pt04 along a high-symmetry path that connects 
the eight time-reversal invariant momenta (TRIMs). The 
majority of the states shown are from Pt and O orbitals, 
whereas the Ca states only appear at above 4 eV. DFT 
predicts the system to be a metal with one single band 
crossing the Fermi level (Ep = 0) being of Pt d x 2_ y 2 
character. 

In addition, we observe the striking presence of a topo¬ 
logical gap at around -4.5 eV. As shown in the light- 
yellow region in Fig. 1, there is a well-established energy 
gap of size 100 meV, the unique feature of which is that 
it is due to two band inversions at T and Z. They are 
characterized by the interchange of the band characters 
(shown as the switch of the green and red colors occur- 
ing at the T and Z points of Fig. 1 from the band below 
to that above the gap in the yellow region.) Specifically, 



FIG. 2. (Color online) (a) The topological surface states show 
two Dirac cones at the (001) surface and no Dirac cone at the 
(111) surface characterizing the weak topological nature of 
Ca 2 Pt 04 - (b) The number of Dirac cones for a given surface 
can be determined from the relative sign of the projected pairs 
of the bulk TRIMs. 

both these inversions involve 0-p z (red circle) and Pt- 
d xy (green circle) bands. This topological gap is purely 
induced by the strong SOC of Pt. In the Supplemen¬ 
tal Material, we calculate the 4 topological Z2 indices 
z/o; (zq, z^2, v 3) as 0; (111), which confirms Ca2Pt04 to be 
a WTI. 

A characteristic feature of a weak TI, as compared to 
the strong TI, is that the topological surface states (TSS) 
do not appear at every surface. This behavior is demon¬ 
strated for Ca2Pt04 in Fig. 2. Two surfaces, i.e. (001) 
and (111), are considered and shown in red and green 
colors, respectively. Evidently, two Dirac cones appear 
at the (001) surface while there is no TSS found at the 
(111) surface. The analysis of the TSS is slightly dif¬ 
ferent from the computation of the four Z2 topological 
invariants and is presented, together with details of the 
slab calculations in the Supplemental Material. 

The predicted topological gap lies deep in the valence 
band. Thus, it will be a challenge in experiments to 
bring the Fermi level down to the topological gap re¬ 
gion. As discussed before, the Ca states basically do 
not contribute to the band inversion, thus, it is possible 
to substitute Ca with other elements with less valence 
electrons to effectively hole-dope the system. Potassium 
substitution may not be sufficient to fully bring the Fermi 
level down to -4.5 eV. However, a subsequent electric gat- 
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ing and/or a direct angle-resolved photoemission spec¬ 
troscopy (ARPES) detection may then become feasible. 

3D topological model: Our concrete example Ca 2 Pt 04 
seems to provide a rich playground for finding new 3D 
topological insulators. This statement holds provided 
that the obviously required inclusion of strong electronic 
correlations in this 5d-material does not destroy the topo¬ 
logical phases identified in the “non-interacting”, i.e. 
DFT + SOC calculations. This problem will be ana¬ 
lyzed next in a 3D topological toy model that resembles 
our material example in the vicinity of the band inversion 
points. Differently from the strong TI, a weak TI can be 
adiabatically connected to stacked layers of 2D TIs [34]. 
This feature motivates us to construct, in 3D, a half-filled 
weak TI model by coupling layers of 2D topological sys¬ 
tems, i.e. 2D BHZ models, with a nearest-neighbor hop¬ 
ping t z along the stacking (z) direction. Subsequently, 
we study the interaction-driven topological phase transi¬ 
tions of this model. 

The “non-interacting” Hamiltonian on a cubic lattice 
is hence given by 

H (k) = [M — (cos k x + cos k y + t z cos k z )]a z <S> I 

+A sin k x • a x (g) S z + A sin k y • a y (g) I , (1) 

where M is the mass parameter, A denotes the SOC 
strength. 0 < t z < 1 is the hopping amplitude along 
the stacked direction (z) perpendicular to the in-plane 
(xy) hopping. When t z is taken as zero, Eq. (1) reduces 
to the 2D BHZ model. This model is different from those 
traditionally introduced for the 3D strong TIs [35-37] in 
that it is still block-diagonal in spin. Note that the Pauli 
matrices cr x ,y,z denote an orbital space degree of freedom 
and the Pauli matrices S x,y,z refers to the physical spin. 

Instead of restricting our discussion to any specific 
choice of the tuning parameters (which will be done in 
our comparison with Ca 2 PtC >4 in the Supplemental Ma¬ 
terial), we study first the Hamiltonian in Eq. (1) from a 
more general perspective and explore all possible phases 
that a topological system, described by Eq. (1), can dis¬ 
play. In the second step, we then want to understand 
how the topological transitions between these phases are 
influenced by electronic correlations. 

The Hamiltonian (1) is time-reversal invariant and at 
the time-reversal invariant momenta (TRIM) k = T^, it 
is simplified to 

H(k = Ti) = [M —(cos k x +cos k y -\-t z cos k z )\cr z <g>I . (2) 

Thus, H commutes with the parity operator P = cr z ® I , 
implying that the eigenstates of Eq. (2) are also the eigen¬ 
states of the parity operator, from which the Z 2 invari¬ 
ants can be easily calculated [11] (see also the Supplemen¬ 
tal Material). As a function of the mass parameter M 
(for a fixed value of interlayer coupling t z ), we find four 
topologically distinct phases, characterized by a different 
number of band inversions. These phases are shown in 
Fig. 3. We identify a weak TI phase, two DSM phases, as 
well as a trivial band-insulating (BI) phase. Notably, also 


the gapless DSM phases are topologically nontrivial [38]. 
Interestingly, in addition to the linear bulk band cross¬ 
ing there are also band inversions at T or at T, X, Z. In 
these phases, the definition of the 3D topological invari¬ 
ants that we used to characterize the weak TI phase is not 
valid any more (because of the gapless bulk). However, a 
2D topological invariant can be calculated for any fixed 
k z plane, as long as a gap remains open in the 2D case. 
Additionally, a mirror Chern number can be defined to 
characterize the nontrivial topology of these phases [38] . 
More specifically, the 3D DSM has its topological origin 
in the 2D Z 2 invariants (for a fixed value of k z ) of the 
corresponding 2D BHZ model. Two 2D topological in¬ 
variants can then be defined at a k z > k c z and k z < k z , 
where at k c z the bulk bands cross. These two 2D topo¬ 
logical invariants differ which is the reason why the bulk 
gap has to close at k z . 
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FIG. 3. (Color online) Depending on the strength of M, the 
3D topological model in Eq. (1) displays four topologically 
distinct phases, characterized by a different number of band- 
inversions as well as the presence (WTI and BI) or the absence 
(DSM) of a bulk gap. 

For the ” non-interacting” 3D layered BHZ model, the 
topological phase transition is driven by the value of M. 
In what follows, we want to show that the same transi¬ 
tions can be triggered by electronic correlations as well. 
To this aim, we add the most general local interaction 
term to Eq. (1), i.e. 

Hjj U ^ ^ ^ ^ ^ ^ia<j^i/3a 

i,a i,a<(3 

T(£7 T) ^ ^ ^iaa^i/3a 
i,a</3,<j 

( C lo;t C */3t C W C ^4 + C la^ C lai C i/3i C i/3t^ (^) 

i,ay£(3 

Here, i is the site index of a 3D cubic lattice, <a, (3 span the 
orbital basis of this Hamiltonian. A rotational invariant 
form of Hu is chosen by us, i.e. U' = U—2J and J is fixed 
to a rather large but not untypical value for transition- 
metals, i.e. U/ 4. 

Evidently, Eqs. (1) and (3) comprise the interacting 
many-body Hamiltonian for a TI that can demonstrate 
various topological phases as shown in Fig. 3 for U = 0. 
In addition to this non-interacting case, there is another 
limit where an explicit solution of the full Hamiltonian is 
known. When U 00 , the paramagnetic ground state 
will be a Mott insulator, which destroys all nontrivial 
topological phases present at weak interactions. 

Before discussing the results in Fig. 4, we briefly de¬ 
scribe the influence of electronic correlations on TIs 
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FIG. 4. (Color online) The interacting phase diagram of the 
3D topological model based on Eqs. (1) and (3). Five differ¬ 
ent phases are identified. At sufficiently large M and U, the 
conventional band insulator and the Mott insulator phases 
dominate. In between, the competition between Coulomb in¬ 
teraction and topological order gives rise to interesting phases 
such as a region of a weak TI phase and two regions of DSM 
phases. The topological surface states (TSS) with parameters 
given by the corresponding colored symbols can be found in 
the Supplemental Material. As we carefully explain there, the 
phase boundaries are subject to an approximation for large 
values of U . 



and explain how the transition boundaries in Fig. 4 
have been obtained. Similarly to their non-interacting 
counterparts, the topology of correlated TIs can be de¬ 
fined by means of an “effective Hamiltonian” h e ^{k) = 
— G~ 1 (k^uj 0) [39-41], which, in the DMFT language, 
reads h e s(k) = H( k) + E(o; -A 0). For any nonzero value 
of M, the imaginary part of the self-energy always extrap¬ 
olates to zero for oj 0. The influence of correlations 
on the electronic structure is therefore mainly encoded in 
5ftF(co> —» 0). Thus, the transitions between the BI, the 
DSM and the WTI phases in Fig. 4 are determined simply 
from an inspection of the topology of i7(k) + — >> 0). 

In the Mott phase, we have however used the even more 
clearcut condition n a = np = 1 [22]. More details about 
the accuracy of the transition lines in the larg e-U part 
of the phase diagram can be found in the Supplemental 
Material. 

The first important observation about Fig. 4 is that 
the topological nontrivial phases extend to a finite value 
of U. Thus, the topology of such a 3D topological sys¬ 


tem is stable against a “certain amount” of electronic 
correlations, which may pave the way for the experimen¬ 
tal realization of correlated TIs. With the increase of 
U , each topologically nontrivial phase extends towards 
a larger value of M. This behavior can be understood 
because the Coulomb interaction modifies the mass pa¬ 
rameter M via 5RF(co> —0). In particular, as an effect 
of the nonzero value chosen for the Hund’s coupling J, a 
larger value of M is required to maintain the phase. 

The second important observation is that, for a fixed 
value of M, the Coulomb interaction U may generate 
transitions between different topological phases. This is 
again due to the inter-orbital interaction which favors 
single orbital occupations. The colored spade, diamond 
and heart symbols in Fig. 4 correspond to three topolog¬ 
ically distinct phases with the same M, but at different 
U values. Among the phase transitions driven by the in¬ 
teraction £/, the first one is of particular interest, i.e. it 
is the transition from a trivial band insulator to a DSM 
and to a weak TI, see Fig. 8(a), (b) and (c). This ob¬ 
servation implies that the Coulomb interaction does not 
always destroy but sometimes also induces topological 
phases. 

In summary, we have investigated the stability of 3D 
topological phases in the presence of electron-electron 
interactions. Our study is motivated by the discovery 
that, on the basis of a DFT band structure calculation, 
the hole-doped compound Ca 2 Pt 04 is a weak TI with 
the band inversions induced by the Pt -d xy and O —p z or¬ 
bitals. This finding in principle opens up a new direction 
in the search for 3D interacting TIs, provided that the 
expected local correlations in this 5 d transition-metal do 
not destabilize the topological nature. To verify this, we 
have modelled the DFT topological band structure as 
the underlying “non-interacting” band structure of a 3D 
generalization of the BHZ Hamiltonian, which has then 
been augmented by the most general on-site two-orbital 
electron-electron interaction. For small to intermediate 
interactions, topological phases as well as transitions be¬ 
tween them are induced by the interactions, in particular, 
a novel DSM to weak TI transition. For very strong in¬ 
teractions, the system enters the expected “trivial” Mott- 
insulating phase. 
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Supplemental Material 


I. CA2PTO4, A HOLE-DOPED WEAK TI 


As a concrete example, we consider the K 2 NiF 4 -type structure of Ca 2 Pt 04 with space group I4/mmm. The lattice 
parameters taken from an ab-initio geometry optimization [42], are a = 3.87 A and c = 11.42 A. The density-functional 
theory (DFT) calculations were carried out within the full-potential linearized augmented plane-wave (FP-LAPW) 
method [43], implemented in the package WIEN2k [44]. K max RMT = 9.0 and a 10 x 10 x 10 k-mesh were used for the 
ground-state calculations. R mt represents the smallest muffin-tin radius and K maa; is the maximum size of reciprocal- 
lattice vectors. The spin-orbit coupling is included by a second variational procedure. The generalized gradient 
approximation (GGA) potential [45, 46] is used in all calculations.. The surface electronic structures are further 
calculated using the maximally localized Wannier functions (MLWFs) [47], employing the WIEN2WANNIER [48] 
interface. The MLWFs are constructed in a non-self-consistent calculation with an 8 x 8 x 8 k-mesh. 


(a) band-inversion 



(b) Weak TI 0;(111) 

r + 



FIG. 5. (Color online) (a). The presence of spin-orbital coupling induces a band-inversion between O -p z (red line) and Pt -d xy 
(green line) bands, which characterises the nontrivial topology of this system, (b) A topological invariant can be defined for 
the topological gap shown in the light-yellow region of Fig. 1. It is found, that this gap possesses a weak-topological nature. 
The +/- sign associated with each high-symmetry point is the parity products for the bands below the topological gap at this 
TRIM. In the coordinate system defined in Fig. 2, the topological invariant is given as 0; (111). 


As shown in Fig. 5(a), without SOC,the O -p z (red line) and Pt -d xy (green line) bands overlap at T and Z. In 
the presence of SOC, the hybridization of these two bands opens an energy gap, which leads to the inversion of the 
two bands at those two TRIMs. The appearance of an even number of band-inversions in a 3D system can give 
rise to a weak TI [34]. Thus, Ca 2 Pt 04 can become a weak TI given that the Fermi level is shifted down into the 
topological gap around -4.5 eV. To more accurately verify the topological nature of this material system at the DFT 
level, we calculate the four topological invariants z/ 0 ; (^,^ 2 ^ 3 ) through the examination of the parity products (Si) 
of the occupied states at all eight TRIMs (see Eqs. (4), (5), (6)). As there is a well-defined energy gap for the band- 
inversions, the separation of the occupied and unoccupied bands at this energy level certainly justifies the definition 
of the topological invariant. Ca 2 Pt 04 respects inversion symmetry. Thus, we can simply identify the parities of all 
occupied bands below this topological gap to calculate the topological invariant. The products of the parities for all 
occupied bands below the topological gap ( i.e . Si in Eq. (5))are marked as +/- sign in Fig. 5(b) for each TRIM. The 
coordinate system of the BZ is chosen to be the same as that in Fig. 2(b). The four Z 2 indices z/q; (zq, zq, zq) are then 
easily calculated as 0; (111). The strong topological index zq = 0, while the other three weak topological indices are 
nonzero which, thus, confirms Ca 2 Pt 04 to be a weak TI. 

Next, we summarize the topological surface states (TSS) calculation of Ca 2 Pt 04 , as displayed in Fig. 2. Taking 
the (001) surface as an example, the surface Brillouin zone (SBZ) has four TRIMs A a , i.e. 1T, 2A, 1M. They are 
the projections of pairs of bulk momenta T a i, T a 2 , that differ by G c /2, onto the (001) plane. As a result, the relative 
sign at the SBZ TRIMs ir a are given as the product of the corresponding M a , i.e. ir a = S a iS a 2 . For example, T is 
the projection of the bulk TRIMs T and A on the (001) surface (see Fig. 2(b)), 7Tp is then computed as £pAv = —1. 
Similarly, one can derive the relative signs at = — 1 and 1 = 1. Thus, there will be two Dirac cones at T and M, 
while there is no TSS at X. In Fig. 2(a), the TSS are displayed from a slab calculation of 8 nm thickness in the given 
direction. In practice, this is achieved by projecting the DFT Bloch bands to the Pt-d and O-p orbitals, which gives 
rise to a tight-binding model with only the target orbitals as a basis. Thus, compared to the full DFT parameter 
space, this basis set is much smaller and, thus, it becomes feasible for a larger slab calculation. In a similar fashion, at 
the (111) surface, one can understand that the momenta 2 N project onto A, 2A onto M, T and Z onto T. Thus, at 



all four surface TRIMs, the relative signs ir a are positive. Consequently, there is no TSS at this surface, see Fig. 2(a) 
for comparison. 


II. A 3 D TOPOLOGICAL MODEL 


By tuning the parameters M, t z and A in Eq. (1), a qualitatively similar electronic structure with a weak topological 
insulator phase - like the one of Ca 2 Pt 04 - can be obtained. In Fig. 6 (a), the bulk electronic structure displays two 
band inversions at Y (0,0,0) and Z (0,0,7r) with M = 1.5, t z = 0.2, and A = 0.3. These two inversions, when projected 
onto different surfaces, can give rise to different numbers of TSS, i.e. a similar situation like what we observed in 
Ca 2 Pt 04 (see Fig. 2). In Fig. 6 (c) and (d), (100) and (001) surfaces are taken as an illustration. At (100), Y and 
Z project into Y and X separately, thus, there are two Dirac cones in the surface Brillouin Zone (SBZ) of the (100) 
surface. In contrast, at the (001) surface, there is no TSS, as Y and Z project to the same Y. 




FIG. 6. (Color online) (a) Bulk electronic structure of the weak TI model defined in Eq.(l) along a high-symmetry path of a 
cubic lattice with the BZ plotted is in (b). Two band-inversions (presented as the interchange of color) at Y and Z give rise to 
topologically distinct surfaces: i.e. at the (100) surface (c), there are two Dirac cones at SBZ momentum points Y and X , at 
the (001) surface (d) there is no Dirac cone. 


In 3D, there are eight TRIMs. This leads to four independent Z 2 invariants for a gaped system. One of these 
invariants, z/q, can be expressed as the product over all eight points, 
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(-ir = IR * ( 4 ) 

i=1 

where Si is the product of the parity eigenvalue for all occupied states at a TRIM In terms of the Hamiltonian (2) 
which has only one degenerate valence band, it take the following simple form [11]: 

Si = -sgn(tf (k = r0) . (5) 

The other three invariants are given by products of four <Vs, for which Yi reside in the same plane, 

(-ir = II S i=(n in2 n 3 ) ■ 

nfc=l;n J *^fc=0,l 


( 6 ) 
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The Hamiltonian (2) at these momenta reads: 


H(T) = [M - (2 + t z )]a z 0 I ; 

(7a) 

H(2X) = [M — t z ]a z 0 I ; 

(7b) 

H(M ) = [M + 2 — t z ]a z 0 I ; 

(7c) 

H(Z) = [M - 2 + t z ]a z 0 / ; 

(7d) 

H(2A) = [M + t z \cr z 0 I ; 

(7e) 

= \M + 2 + t z ]cr z 0 I. 

(7f) 


At each TRIM, the Hamiltonian is proportional to the parity operator P = cr z <g> I, thus, they obviously have the 
same eigenstates. Furthermore, at each TRIM, two of these states form one occupied Kramers pair and the two other 
ones correspond to an empty Kramers pair. If M > 2 + t z , at all eight TRIMs, the occupied Kramers pair has an 
eigenvalue of —1. Therefore Si = 1 at all eight TRIMs and the topological invariants are simply 0; (000), this is the 
trivial band insulator (BI) phase shown in Fig. 3. For 2 — t z <M<2 + t z there is a bulk-band linear crossing between 
Z and T (see the left plot in Fig. 7) and the system is a Dirac semimetal. According to the classification for 3D 
Dirac semimetals, the DSM in our model belongs to the first class that is created via a band inversion [38], which is 
topologically nontrivial (see also the topological characterization of this phase in the main text). The appearance of 


M = 1.8, A = 0.3, t z = 0.5 M = 0.2, A = 0.3, t z = 0.5 




FIG. 7. (Color online) The bulk band dispersion for the two DSM phases, (left) When t z < M < (2 — t z ), there is a bulk-band 
linear crossing between T and Z; (right) when M <t z , the linear crossing is between A and X. 

the bulk-band linear crossing is one striking feature of this 3D topological model, which is resulting from the different 
inversions at Z and T. The k z = 0 and k z = tt planes are different only at T and Z, as one is inverted and the 
other one is not. Thus, when changing k z = 0 to k z = 7r, the topology has to be changed from a nontrivial one to 
a trivial one. One way of changing the topology is to close the bulk gap. Since there is a band inversion at T, this 
DSM phase is still topologically nontrivial [38]. The difference to a strong TI is the absence of a finite bulk gap in 
the DSM. Thus, the corresponding topological surface states in the DSM phase are embedded in the bulk bands. If 
t z < M < (2 — t z \ a weak TI phase appears, in which the occupied Kramers pair at both T and Z have eigenvalues 
1. The corresponding topological invariants are then 0; (001). For M < t z a DSM phase appears again, in which the 
bands at the two X points are also inverted. Similarly to the difference between T and Z for t z < M < (2 — t z ), now, 
the topology of the k z — 0 and k z = tt planes only differ at X and A. Thus, by changing k z from 0 to 7r, we have to 
go through a bulk gap closing between A and X (see the right plot in Fig. 7). The complete phase diagram of the 
Hamiltonian (1) for M > 0 is displayed in Fig. 3. It is symmetric with respect to positive and negative values of M. 
For negative values of M, bands at more TRIMs will be inverted. However, the different phases remain. Thus, in the 
following many-body investigations, we only focus on the positive choice of M. 

When U ^ 0, we have to solve the Hamiltonian at finite U numerically. Here, we employ dynamical mean-field 
theory (DMFT) [49] with the continuous-time hybridization-expansion (CT-HYB) [50, 51] quantum Monte Carlo as 
an impurity solver. The CT-HYB solver is based on the implementation discussed in previous works of us [52, 53]. To 
illustrate the correlation effect, we set the SOC strength to A = 0.3 and the inter-layer coupling to t z = 0.5 throughout 
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the calculations and explore the interacting phase diagram as functions of M and U. The phase boundaries in Fig. 4 
are determined from the effective Hamiltonian H( k) + 5RE(u; —>■ 0), as explained in the main text. The transition lines 
to the DSM phases are more difficult to be located, as they involve gap closings, which break down the definition of 
the 3D topological invariants. A full analysis would actually require a study of the local spectral function obtained 
from both the real and the imaginary parts of E(cj) or a careful examination of the temperature scaling of the 
Green’s function. Here we have taken a simplified route and looked only at the gap closing induced by the rigid 
shift —)> 0), which is accurate mostly in the weak -U region of the phase diagram. We note that the precise 

location of the phase boundaries is not of our prime interest here, as we are interested in the general trend of the 
phase transformation as a function of U. 

In Fig. 8, the topological surface states (TSS) of four topologically distinct phases are displayed with parameters 
given by the corresponding colored suit symbols shown in Fig. 4. The TSS are calculated for the (100) surface (left 
plot) and (001) surface (right plot) in a 40-layer slab. For simplicity and to have a better low-energy resolution, we 
have used the effective one-particle Hamiltonian h e s(k) with only the real-part of the self-energy SRE(u; —>> 0) in the 
calculations. At uj ^ 0, the nonzero value of 9E(u;) broadens the electronic bands and smears out the details of the 
TSS at large frequencies. However, the topology of the system is fully determined by h e ff(k, uj = 0) where ^E(c<;) is 
zero in all four phases displayed in Fig. 8. In Fig. 8(a)-(c), we show the TSS with fixed M — 3 and three different 
values of U. With the increase of interaction, a topological nontrivial phase can be transformed to a DSM by inverting 
the bands at the T point and further to a WTI with invariants 0; (001) by also inverting the bands at the Z point. 
In Fig. 8(d), we show another type of DSM, the corresponding parameters are taken as M = 1, U = 4. Its TSS also 
contains two Dirac cones, but their locations in the SBZ are different from that in Fig. 8(c) (see also Fig. 3 for the 
difference of band inversions). The left four plots in Fig. 8 are the TSS at the (100) surface, which all contain a 
well-defined band gap, however this is not the case for the (001) surface. At the (001) surface, the two DSM phases 
are gapless. In combination with the left plots, one can easily understand that the TSS in the two DSM phases are 
embedded in the bulk states. 



FIG. 8. (Color online) (left) The topological states at the (100) surface with (a) U = 4.0, M = 3.0; (b) U = 7.0, M = 3.0; (c) 
U = 8.0, M = 3.0; (d) U ~ 4.0, M — 1.0. From (a) to (c), the topological phase transitions are driven by the increase of U. 
(right) Same as the left plot, but for (001) surface. At the two DSM phases, the bulk gaps are closed, thus the TSS at these 
two phases are embedded in the bulk bands. 



















































